NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 




TECHNICAL MEMORANDUM 1260 


EXACT SOLUTIONS OF EQUATIONS OF GAS DYNAMICS 

By I. A. Klebel 

Translation 

“Primer Tochnogo Reshenia Ploskoi Vikhrevoi Zadachi Gazovoi 
Dinamiki.” Prikladnaya Matematika i Mekhanica. Vol. XI, 1947. 


^NACA^ 

Washington 
June 1950 


AFMOC 

Tirr ■ 1 j • n ^ s n \j 

I ‘ . < innl 


2 If ff/ 


- 1 1 


TECH LIBRARY KAFB 



TECH LIBRARY KAFB, NM 


NATIONAL ADVISORY COMMITTEE FOE AERONAUTIC 


Mill 

OlMMbTS 


TECHNICAL MEMORANDUM 1260 


EXACT SOLUTIONS OF EQUATIONS OF GAS DYNAMICS* 

By I. A. Kiebel 

The equations of the two -dimensional stationary problem of gas 
dynamics are of the form 

Su Su 1 Sp 

ox ay p Ox 


Sv St 1 Sp 
U Sx + V Sy = “ p Sy 


^ + *EL „ o 
Sx by 




( 1 ) 


J 


where u «nfl t are the components of the velocity along the 
coordinate axes x and y, respectively, p is the pressure, 
p the density, and K= c p /c T , the ratio of specific heats. 

The equation of continuity permits the introduction of the stream 
function from the equations 


Pu = 


S\J/ 

dy 




pv = - — 


Sx 


( 2 ) 


J 


The last of equations (l) gives 


1 

K 


P* 0l0 


(3) 


*"Primer Tochnogo Reshenia Ploskoi Vikhrevoi Zadachi Gazovoi 
Dinamiki." Prikladnaya Matematika i Mekhanica. Vol. XI, 1947, 
pp. 193-198. 
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■where d is a certain function only of The first two of 

equations (l) give the Bernoulli law 


u 2 +v 2 


w-1 

K 


+ ~ - io (*> 

IV X 


(4) 


where i Q 
is 


is a function of and the equation for the vorticity 


br 5 u 
bx by 



K-l \ 
K K d* ] 

K-l p d^/ 


(5) 


For the solution of the vortex problem in which at least one 
of the derivatives di Q /d^ and d$/d ^ is different from zero, it 
is convenient to pass from the variables x and y to the vari- 
ables x* = x and 'l/. Equations (2) and (5) then assume the form 
(See, for example, reference 1.) 

by v 

dx* * u 

<k « JL 

B'J/ pu 

5v 
3x* 

The problem reduces to the determination of the five functions 
u, v, p, p, and y of x* and from equations (3), (4), 
and (6) . 1 

By using the last of equations (6), the function X i 0 intro- 
duced with the aid of the equations 


. b \ |/ 


> ( 6 ) 
J 
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The first two equations of equation (6), on the "basis of equa- 
tions (3) and (4), then assume the forms 

K-l . _ 1 


Sy 

ax 


ax 

bil 


2io - 


2K_ 

K-l 


4 IE 




( 8 ) 


a* # 


i 

K 


2i 2K_ fax’ 

21 *~ 4 vs; 


K-l 

K. 


\ \2 

ax 

vow 


( 9 ) 


where the asterisk on x has been dropped. 

By differentiating equation (8) with respect to and equa- 
tion (9) with respect to x, y can be eliminated and a single 
equation for the function X obtained. 


In order to obtain an example of the exact solution of the 
system of equations, i Q is set equal to a constant: 

K-l 

X * - H 010 x K+1 (10) 

where H is a certain function of to be determined. Equa- 
tions (8) and (9) then give 



( 12 ) 
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When equation (11) is differentiated with respect to v|/ and 
equation (12) is differentiated with respect to x, terms with 
the sa me degree of x are collected and after simple trans- 
formations there is obtained 

_1 

* (*=1 H ) " . h " 

\j«+l J 2 

H* (H ,2 +KH H")* = - (K-l) H" (H ,2 +KH E")J 

Successive integration qf the second of equations (13) gives 


/ (13) 



From the last equation, the relation between ^ and H can 
be found with the aid of quadratures. It is more convenient, how- 
ever, to introduce a numbering of the streamlines directly with 
the aid of H' and not 


Thus, equations (14) give 1 as a function of H’ and, using 
equations (13), permit finding $ in terms of H' in the form 


t = 


K-l 
2 K 


02 


K-V 

K+l 


1-K 

K 




2K 

(C+l 


(15) 


Substituting in equation (11) and integrating yields, for the 
streamlines. 


y = 


U 


H'x 


K-l 

k+t 


h - 


-2 


K-l 1 


,1-K, * +1 l 2 


CjiH- ‘X 


dx + F (H* ) 


(16) 
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Comparison with, equation (12) 'shows that F = constant and 
without loss of generality can "be taken as F = 0. 

i 

For the determination of u, v, p, and p the following 
relations are written: 


u = 


2i 0 - C 1 H ,1 " lt x 


1-K 

1+K 




H'x 


K-l 

K+l 


( 17 ) 


P 


C 2 



K 2k 
1+ K 1+ K 
X 


p = 


2K 

K-l 



2 K+l 
K+l 


E* 



2 

1+K 


( 18 ) 


( 19 ) 


For simplification it^is convenient to introduce the nondimen- 
sional variables x and y and in place of H' introduce t| 
with the aid of the following equations : 


x = Lx 
7 = 1-7 




H 


1 1+K 


'1 ~2 j 


( 20 ) 


where 


_ 1 K+l _1_ 

L = (2i 0 ) 2 K_1 C^" 1 (21) 


Equations (16) and (17) then assume the following simple form: 
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x = Tjm (y) 


2 oj /i K +1\ 

u = 2i D (1-m ) 


K-l 

1 "k+1 


Y = V 2i o n m = V 2i o | m 


1 Kfl 


gives 


If K » 1.4, the integral in equations (22) is evaluated and 

i 

y « 6 - 1 1 + | ( ^£-1) + i ( ^-l) 2 (24) 


and the streamlines can therefore he easily constructed. 

Figure 1 shows the streamlines for t) = 1, 2, 3, . . . , 19. 

This motion possesses both supersonic and subsonic velocities, 
for the line of transition (shown dotted in fig. l) is obtained If 


u2 * t2 - 2 OT 


_ g ±L -2 — 

, K+1 1 K+1 K-l 

1-m + -TS" m = — — 

T ) 2 * +1 


o K-l 

- K+1 - K # f 1 ' ^ 


► 
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It is possible, without difficulty, to construct the character- 
istics in the x,y- and u,v -planes . Instead of the equation of 
the epicycloid in the u,v-plane, in the case considered a more 
complicated equation arises . 

In the vorticity problem, along the characteristic there 
occurs 


dp =F 


ctg a 
v 


dw 


± 


sin a cos a 

ic^L 


d log £ 


(26) 


where a is the Mach angle, w the magnitude of the velocity, and 
p the angle of inclination of the velocity to the x-axiB . Equa- 
tions (23) yield 


*1 


2 


2i c -u 


2i Q -v 2 cos 2 p 
w 2 sin 2 0 


(27) 


But 

equation 


depends only on H ' , that is , on 
(15), the following is obtained: 



According to 




K -1 
2k 



2k 

(tj 2 -1)K+1 


K-l 

2k 


1-K 

n *zk\ K 
°2 STl) 


2i -w cos p 


K-l 2k 

(v 2 sin 2 p) K+l (2i 0 -w 2 )*+l 


(28) 


Substituting this value of £ in the right side of equa- 
tion (26) yields the differential equation for the characteristics 
in the p ,v-plane . 

Finally, the question arises whether it is possible to obtain 
such vortex motion by transition through a surface of strong dis- 
continuity of some other kind but with irrotational motion. This 
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problem may "be answered in the affirmative. For on a surface of 
discontinuity there must, among others, he satisfied the relations 


P+ 

P_ 


K-l 2 K 
K+l K +1 


p+ e + 



*p+ \ p + * + 

1 T 1 


p + 0, 


p+ 




(29) 


where p + , p + , and 9 + are the pressure, density, and velocity 
of propagation of the surface of discontinuity on one side of this 
surface and p_ and p_ are the pressure and density on the other 
side. The magnitudes p + and p + may he taken from the vortex 
motion. The magnitude is found in terms of the elements des- 

cribing the motion and in terms of the inclination of the surface 
of discontinuity. Finally, p_ and p_ can he connected hy the 

relation p_ ' = where is a constant (up to a oertain 

degree of arbitrariness ) magnitude. (At the left of the surface of 
discontinuity the motion is irrotational . ) 


\ 


Inasmuch as 


then 


where 



2K 

tt+1 








(30) 


(31) 


( 32 ) 


-9 = V cos S + u sin & 

+ 4 * + 


Inasmuch as 
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where 5 is the angle between the normal to the surface of dis- 
continuity and the x-axis (the normal is directed toward the 
"positive" region) M may he expressed in terms of known magni- 
tudes and 8 . 

By expressing tan 8 in terms of the derivative of x 
with respect to rj along the surface of discontinuity and using 
equations (22), (18), and (19), the differential equation for 
determining the surface of discontinuity (K = 1.4) is obtained 
after simple transformations : 



equation (31) in which the constant G- is, to a great extent, arbi- 
trary. The velocity in the "negative" region will be determined 
first on the surface of discontinuity and then extended on the 
negative region by the usual graphical method of Busemann. A model 
of the motion about a contour having an angle is obtained. 

The surface of discontinuity extends out from the angle and 
on passing through the angle the motion reconverts to the rotational 
motion herein considered. 

A solution analogous to that developed can also be obtained 
for the problem with axial symmetry. Taking for the independent 
variable the distance r* = r from the axis of symmetry z and 
the stream function \J/ yields the relations 

Sz Y z 
^ ~ v r 

Sz 1 
Sr* - r*pv r 

where v r and v z are the velocity components. As before, 

1 

p K = & (^)p 


~\ 

> (34) 

J 
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Bernoulli's lav vill have the form 


1 

2 



+ 


K-l 


K-l 

IT 

V 


i 


o 


The equations of Euler give 


dv z dr*p 
dr* “ ~W 


The function X(r*,^) can therefore he introduced from the 
conditions 


T z - - 


5? 


p = 


_ 1 c>X 

r* dr* 


(35) 


Equations (34) and (35) nov permit writing (the asterisks on 
are dropped) 


dz _ dX JL_ 
d'l' = dty v r 


K+l 

h.r * m 

dr 




(36) 


where 


k-l 




By eliminating z from equations (36), a single equation for 
the determination of the function X is obtained. Particular 



NACA E4 1260 


solutions, analogous to the solutions in. the first case, can he 
constructed "by seeking X in the form 


_ £ 

X « - r K+1 H(i|r) 


Translated "by S. Reiss 
National Advisory Committee 
for Aeronautics. 
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Figure 1. 
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